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Abstract 


This progress report presents the results of an investigation focused on parameter identification for the 
NASA F/A-18 HARV. This aircraft was used in the high alpha research program at the NASA Dryden Flight 
Research Center. In this study the longitudinal and lateral-directional stability derivatives are estimated from flight 
data using the Maximum Likelihood method coupled with a Newton-Raphson minimization technique. The objective 
is to estimate an aerodynamic model describing the aircraft dynamics over a range of angle of attack from 5° to 
60°. The mathematical model is built using the traditional static and dynamic derivative buildup. Flight data used 
in this analysis were from a variety of maneuvers. The longitudinal maneuvers included large amplitude multiple 
doublets, optimal inputs, frequency sweeps and pilot pitch stick inputs. The lateral-directional maneuvers consisted 
of large amplitude multiple doublets, optimal inputs and pilot stick and rudder inputs. The parameter estimation 
code pEst, developed at NASA Dryden, was used in this investigation. Results of the estimation process from 
a =5° to a =60° are presented and discussed. 
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Symbols 


a,, = normal acceleration, g 

a x = longitudinal acceleration, g 

ay = lateral acceleration, g 

b = wing span, ft 

c" = mean aerodynamic chord, ft 

= aerodynamic coefficient where i = (N,m,A), rad 1 or deg 1 
E[ ] = expected value 

g = gravity acceleration, ft/sec 2 

I = moment of inertia, slug-ft 2 

m = aircraft mass, slugs 

p = roll rate, deg/sec 

q = pitch rate, deg/sec 

"q = dynamic pressure, lbs/ft 2 

R = conversion from radians to degrees (57.3) 
r = yaw rate, deg/sec 

S = wing planform area, ft 2 

V = velocity, ft/sec 

W/T - wind tunnel 

x = sensor location, ft 

y = computed aircraft response vector 

z = measured aircraft response vector 

Greek 

a = angle-of-attack, deg 

(3 = sideslip angle, deg 

5 = control surface deflection, deg 

0 = pitch attitude, deg 

£ = parameter vector to be estimated 

1 = summation 

4> = roll attitude, deg 

v = gradient 

Subscripts 

A = axial force 

a = aileron 

ba = basic airframe 

D = drag force 

dht = differential horizontal tail 
e = elevator 

L = lift force 

1 = rolling moment 

lef = leading edge flap 

m = pitching moment 

N = normal force 

n = yawing moment 

pv = pitch vane 

r = rudder 

sa = symmetric aileron 

t = true parameter value 

tef = trailing edge flap 

Y = lateral sideforce 

yv = yaw vane 
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Introduction 


The NASA High Alpha Technology Program (HATP) was initiated with the goal of obtaining a better 
understanding of aircraft aerodynamics, control surface effectiveness and airflow phenomena at high alpha flight 
conditions. Successful achievement of these goals are envisioned to allow expansion of flight envelopes, 
maneuverability enhancement and flight safety for future aircraft. The NASA F/A-18 High Alpha Research Vehicle 
(HARV) project is one of the flight test aircraft involved in this program. 

The research work presented in this document is related to parameter identification (PID) for the HARV. 
The goal is to determine an aerodynamic model which accurately describes the dynamics of the HARV over a wide 
range of angle of attack conditions. This process allows the validation of wind tunnel data and results obtained from 
other PID methods. In addition, the development of a more accurate model is of great importance for updating 
flight simulators and onboard flight computers by enabling aircraft designers to better emulate the true system 
dynamics. Results from PID analysis may also be used as a tool for monitoring aircraft stability as well as control 
surface effectiveness while investigating the expansion of the flight envelope into the high alpha regime. 

This report is related to the longitudinal and lateral-directional PID results for the HARV over a range from 
5° to 60° for the angle of attack. The derivatives are extracted from flight data using the Maximum Likelihood 
(ML) method. The ML cost functional is maximized through the minimization of a quadratic cost function which 
contains the differences between measured and computed system responses. This minimization is achieved using 
the Newton-Raphson (NR) technique. The parameter estimation program (pEst), developed at NASA Dryden, is 
implemented in this research. 

Flight data analyzed includes a variety of longitudinal and lateral -directional maneuvers obtained from the 
HARV flight testing activities. Particularly, the longitudinal maneuvers examined include single surface input 
multiple doublets (SSI MD), optimal longitudinal (OLON), 3211 frequency sweeps, and pilot pitch stick steps 
(PPSS). The lateral-directional flight data included single surface input multiple doublets (SSI MD), optimal lateral 
(OLAT) and pilot yaw/roll rudder and stick steps (PYRS) obtained from HARV. 
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Chapterl: Maximum Likelihood Method and Newton Raphson Algorithm 

The ML method coupled with a NR minimization technique has been one of the most successful PID 
approaches for several years. The effectiveness of this approach is well documented and excellent results have been 
achieved for a large variety of aircraft, as presented in Refs. [1 ,2, 4, 5, 7, 8, 9]. It is known that the estimates obtained 
using the ML method have three asymptotic properties: they are unbiased, approach a Gaussian distribution and have 
the lowest possible variance 6 . By asymptotic it implies that the aforementioned estimate properties are true if intinite 
data time is available. However, these estimate properties are best approximated if the data time is long enough 
(ie. a couple of periods of the lowest system natural frequency) 6 . In the case of an aircraft, "long enough” would 
be defined as a few periods of the phugoid mode. 

In general, the objective is to maximize the probability that the computed system responses, based upon 
a set of estimated stability derivatives, are representative of the true system dynamics. This conditional probability, 
denoted as P(z/£), is also known as the likelihood function. In this analysis it is assumed that all state and 
measurement noise may be described as a Gaussian, white sequence with zero mean. This allows the likelihood 
function to take on the following form: 


P(z/i) =- 


/jVM 

2 2 Vr 1 ! 




( 1 . 1 ) 


where the Gaussian probability density function P(z/§) represents the ML cost functional to be maximized. Here 
Wi represents the response error weighting matrix. The number of discrete time points in the maneuver and the 
number of responses in the cost function are denoted by n, and n, respectively. J(£) represents the quadratic cost 
function to be minimized in order to maximize the likelihood function. This cost function is defined as follows: 


= Y ([z(t k )-y(c k )] T H 1 [z(t k )-y(t k )]) 11 ’ ] 

A n t n z JPa 

By expressing J(^) as its second order Taylor Series expansion and setting VJ(£) equal to zero the NR technique 
generates an expression for updating the parameter vector at each iteration in the estimation process: 

5i-i= 5i - [(VV^))]' 1 (1 ‘ 3) 
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where the first gradient of J(Q with respect to £ is given by: 


E[ [V {y (t Jc)] r (z(t k )-y(t k ) )] (1 - 4) 

n z m 

The Gauss-Newton approximation of the second gradient of J(Q with respect to % is as follows: 

vy<W-i ^ w, v t y(t*)] (1 - 5 > 

A block diagram of the overall estimation procedure is shown in Figure 1 . As illustrated, a predetermined 
set of control surface inputs are applied to the true aircraft system and the responses are recorded. The recorded 
control surface inputs and resulting aircraft responses are then read into the pEst code 3 . 

The aircraft mathematical model within pEst contains the polar coordinate form (a,p,V) of the 6 degree 
of freedom (DOF) non-linear aircraft equations of motion with the state equations being integrated using a fourth 
order Runge-Kutta numerical scheme. In addition, when computing estimated responses pEst also accounts for 
sensor locations with respect to the aircraft center of gravity (CG) to ensure that the cost function is comparing the 
same "local" responses recorded by the sensors. Most importantly, pEst is well-suited for user interaction allowing 
the analyst to determine which stability derivatives are to be estimated and which states or responses are to be 
included in the quadratic cost function. The user may also substitute measured flight data in place of those states 
or responses that are not computed within pEst. In addition, the diagonal weighting matrix (W,) is easily set by 
the user to indicate which response errors are to be considered more important in the estimation procedure 

The recorded control inputs are then applied to the mathematical model of the aircraft within pEst and the 
differences between the measured and estimated responses are determined. The NR algorithm then determines the 
updated parameter vector to be used for the next iteration This procedure continues until a desired convergence 
criteria for the change in the quadratic cost function AJi \) from one iteration to the next is met resulting in the final 
ML estimates of the aircraft model parameters. 

The code is also capable of computing the corresponding Cramer-Rao bound for each derivative. This 
bound is a computed minimum standard deviation of the estimates which can be used as a measure of their 
goodness. With the assumption of asymptotically unbiased estimates the Cramer-Rao inequality may be written as: 
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( 1 . 6 ) 


o($) = JeUFTT) TTTJ’TU * 

The assumption of Gaussian, white system and measurement noise allows the Fisher information matrix to 
be approximated as the Hessian of the quadratic cost function J(£): 

W(5 t ) “ -z-V E[ [V t y(t*)] r V { y(t*)] = V* C J<W (1 ' 7) 

n t n z k^L 

This derivation is shown in Refs. [1,6]. In general, the most practical use of the Cramer-Rao bound is as a measure 
of scatter in the final estimates for the given flight data 6 . 
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Chapter 2: Estimation of the Longitudinal Aerodynamic Parameters 


2.1 Mathematical Modeling of the Longitudinal Aircraft Dynamics 

Modeling of the aircraft longitudinal dynamics within pEst begins with the polar coordinate form (a,p,V) 
of the classical 6 DOF, non-linear equations of motion as seen in 4,9,11 . The body axis longitudinal equations of 
motion are shown below: 

& = g-tanp (p cosa+r since) - 

— (cos6cos4>cosce+sin8sintt) - (2.1) 

Vcosp /nVcosp 

= M y + [pr(I 2 -I x ) *(r 2 -p 2 ) I^+qrl^-pqly,] (2.2) 

0 = q cos<t> - r sincfc (2.3) 

V r =-_2?c D +g r (cos<t)cos0sinacosp+sin<t)cos0sinP - 


3 in 0 cosacosp) +— cosacosp 

777 


(2.4) 


where : 

C L -C N cos (a) -C A sin (a) (2.5) 

C D =C A cos (a) +C^sin (a) (2.6) 

In addition, adjustments were built into pEst to include the effects of angular rates and accelerations to the computed 
data to account for accelerometers which are non-coincident with the aircraft CG. The additional estimated 
parameters a^ and a* Bil# are necessary to account for accelerometer instrumentation biases as shown: 
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= |f C w - ^ [ (x a - Xcg ) H ♦ W a -y C9 )P) ~ 


(V 2 + P 2 ) + a. 


nBiae 


(2.7) 


- -&c A - ~^[(z a -z cg )4 * (y..-y«)*] 


mg gi? 


a x 7 C ?' 


S* 


- (x ax -x c? ) (g 2 +r 2 ) + a 


xBiaff 


( 2 . 8 ) 


Similarly, the flow angle measurements for a recorded by the aircraft instrumentation are also affected by angular 
velocities and upwash if the flow vane is not coincident with the aircraft CG. Consequently, the computed flow 
angle a is adjusted within pEst as follows: 


a = k„ [a+ (x.-x^) (y a ~y cg ) 


(2.9) 


Sensor and CG locations above, behind and to the right of the reference point are considered positive 3 . The 
reference point is located approximately 5 ft in front of the aircraft nose. 

The linear buildups used to model the total, non-dimensional aerodynamic force and moment coefficients 

within pEst are given by: 

c w = C No * C N a * C N Je + C Nti J 1 ef* C„ it J tef+C^Jsa * 




( 2 . 10 ) 


c m = C Wo * C m a * c m je + 1 ef * tef * C mt Jsa + C m j-^q 


( 2 . 11 ) 


C A = C Ao ♦ Cj^a * C^b e * C A(i J 1 ef * C^J tef * C At Jsa * 


( c *>Jnik +c 4^ ) q 


( 2 . 12 ) 
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Note that the static buildup for the moment coefficient does not include the derivatives describing pitching moment 
effects due to the pitch vane. This moment derivative is found in the total pitching moment equation where: 

My = qScCJi-6 1 ( ERPM ( I xe ) ) + * [ C v «pv] R (2.13) 

The term ERPM represents engine revolutions per minute while I« denotes engine moment of inertia. The factor 
of 6 is simply a conversion from revolutions per minute to degrees per second. Note that a standard sign convention 
characterizing conventional control surface motion was employed: positive control surface deflection denoted as 
being trailing edge down. The overall pitch vane deflection was defined as positive for downward deflection of the 
exhaust plume. 


2.2 Estimation Procedures 

Four classes of longitudinal PID maneuvers were available for analysis in this effort. The SSI MD 
maneuvers consisted of independent longitudinal control surface doublets generated by OBES for the trailing-edge 
flaps, symmetrically deflected ailerons and the stabilator. The OLON runs included a series of 7.5° stabilator nose- 
up and nose-down pulses applied by OBES. The 321 1 runs involved an OBES generated frequency sweep of the 
stabilator with a 7.5° amplitude of deflection. Such a maneuver consisted of four stabilator steps that alternate in 
sign with the first being three time units long, the second being two units in length and the last two pulses each 
being one time unit in duration. Each of the aforementioned maneuvers were performed with OBES on and the 
TVCS active. The PPSS runs, consisting of pilot pitch stick step inputs that alternate in sign, were divided into two 
subclasses: pilot pitch stick steps with the TVCS active (PPSSTV) and without (PPSS). Representative plots of the 
control surface inputs for the aforementioned maneuver types may be seen in Figure 2 recorded at an alpha of 20° 
In each class of maneuver investigated, the longitudinal inputs were designed to minimize lateral-directional 
dynamics. All maneuvers were initiated under trimmed, straight and level flight conditions with all flight data 
sampled with a frequency of 40 Hz. The test maneuvers were performed with a fairly consistent CG location 
varying over the range of 24.5% - 26.7% c . 
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Concerns regarding the possible instability and overestimation in magnitude of estimated derivatives due 
to high correlations in the flight data induced a need for model simplification. The presence of correlations among 
static and/or dynamic flight data implies that a dependency exists among them often leading to difficulties in 
extracting model parameters that describe their separate effects on the overall system. Such correlations often exist 
among control surface deflections which are dependently deployed with respect to other control surfaces, or other 
static flight data in the case of flight control system alpha scheduling. Analysis of the HARV flight data indicated 
that, under normal operating conditions, the leading-edge and trailing-edge flaps are deployed using a schedule 
imparting high correlations between each other, as well as angle of attack. The presence of this high correlation 
can clearly be seen in the sample time history of Figure 3a. To avoid biases in the derivative estimates all leading- 
edge flap derivatives (C Naef , and C ASef ) were held constant at their relative wind tunnel values dunng the 
estimation procedure. Similarly, all trailing-edge flap derivatives (C N5wf , C [u5(cf and C A5lef ) were held constant at 
wind tunnel data for all maneuvers with the exception of those in which they were independently pulsed by OBES 
as in the case of the SSI MD maneuvers. 

A similar technique was employed while analyzing the PPSSTV flight data in which high correlations 
between stabilator and pitch vane deflections, as can be seen in Figure 3b, were encountered due to a lack of 
independent control surface inputs. As a result, in an attempt to achieve more accurate estimates of the pitch vane 
coefficients, the stabilator derivatives (C N6e , and C^j were held constant at wind tunnel data dunng the 
estimation procedure. All wind tunnel values were based upon angle of attack and Mach for the maneuvers. 

All mass properties (x,,, y cg , z c ,, I x , I y , I 2 , I u and mass) were taken to be unique values for each individual 
maneuver investigated. These values were obtained from flight dau with an average value of each property over 
the appropriate maneuver time being used in the mathematical model found in pEst. Finally, all lateral states and 
responses found in the aircraft equations of motion wrre replaced with measured flight data throughout the 
estimation process, as was the aircraft flight velocity V 

This investigation produced a set of estimates for the following static and dynamic longitudinal derivatives 
for the NASA F/A-18 HARV: C No , C Na , C N6e , C N4 ^. C.,^. C Nflu , C Nq , C,,*,, C m , C^, Cmspv. CmSsa> 

0*,, C Ao , C Aa , C A5c , C A5tef , C A6pv , C ASs „ C a , In addition, the two parameters a nBui and a* Bus were estimated 
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in order to model the measurement biases for the normal and axial linear accelerometers aboard the HARV. A 
complete listing of the estimated aerodynamic coefficients, as well as which derivatives were held constant at wind 
tunnel values, for each class of maneuver may be seen in Table 1. 

In terms of the longitudinal cost function (J), the objective was to minimize the difference between the 
measured and computed values of a, q, 0, and a, as shown below: 




2 n t es 


a ( t k ) -& ( t k ) 
q{ t k ) -£( t k ) 
0 ( t k ) -0 ( t k ) 

3 x (t k ) S x (t k ) 


A 0 0 0 O' 
0 B 0 0 0 
0 0 C 0 0 
0 0 0 D 0 
0 0 0 0 E 


' a ( t k ) -ft ( t k ) 1 
q{t k ) -Q( C k ) 
0(t*) -0 ( C k ) 
a n (t*) S„(t k ) 
3 x (t k ) -S x (C k )[ 


(2.13) 


where the terms a, q, 0, and a x are the responses measured from the aircraft. The corresponding responses 
computed in pEst are denoted by &, q, 0, and S x . The diagonal matrix within the cost function represents the 
user defined weighting matrix (W,). A trial and error approach was used to determine the best set of weightings 
for each class of maneuver investigated. A listing of the A, B, C, D and E elements of the weighting matrix, as 
was used for each class of maneuver, is given in Table 2. 


2.3 Results of the Parameter Estimation Process 

The results of the investigation are in terms of comparisons of actual and calculated time histories as well 
as estimates of the longitudinal derivatives and a comparison with the wind tunnel estimates. Although estimates 
were obtained for the C N , C m and C A derivative subcomponents this paper emphasizes the subcomponents of C N and 
C m only. This is due to the fact that the maneuvers under investigation did not adequately excite axial acceleration 
(a*) or aircraft flight velocity enough to obtain confident estimates of the axial coefficients. 

An analysis of the actual and calculated time histories reveals that, overall, the minimization technique has 
resulted in an accurate representation of the measured HARV responses even at higher angles of attack. A typical 
sample time history match may be seen in Figure 4 showing time history comparisons between measured and 
computed a, q, 0 and a„ responses at 30° alpha for the SSI MD class maneuver. Of course, accurate fits between 
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measured and computed aircraft responses does not alone imply that an accurate aerodynamic model has been 
estimated. 

The final PID results for the main longitudinal static and dynamic parameters are shown in Figures 5-7 
related to an angle of attack range of 5*- 60”. Each estimate in the plots indicate the maneuver type spawning that 
estimate as well as its Cramer-Rao bound graphically represented as a vertical bar. Recall that the Cramer-Rao 
bounds are equivalent to an estimated minimum standard deviation for a given estimate which, ideally, should 
accurately describe the scatter of the estimates at each angle of attack 6 . Some considerations are needed for a better 
understanding of results in terms of the Cramer-Rao bounds. In fact, recall that the determination of the Cramer- 
Rao bounds assumes unbiased estimates with system and measurement noise being described as a Gaussian, white 
sequence. Considering a case where these assumptions are truly valid the time history residuals would also be 
white, random sequences resulting in a set of Cramer-Rao bounds accurately modeling the scatter of the estimates. 
However, such ideal test conditions are only guaranteed if flight data obtained through computer simulation is being 
investigated. As expressed in Ref. [6], analysis of real flight data may result in biased estimates often caused by 
modeling inadequacies or poor system excitation. Furthermore, colored system or measurement noise will also 
affect the computed residuals As a result, in the case of residuals correlated in this manner, optimistically 
underestimated Cramer-Rao bounds are often calculated causing an inadequate model for the scatter of the estimates. 

A variety of techniques to "correct" the Cramer-Rao bounds are discussed in detail in Ref. [6]. One such 
technique involves multiplying the bounds by a factor, often 5-10, to compensate for the effects of the colored 
residuals. A more precise method in Ref. [6] is based on the determination of a unique correction multiplier for 
each Cramer-Rao bound through examination of the power spectral density for the response residuals generated by 
the estimates. Note that the Cramer-Rao bounds presented herein are not adjusted. 

In order to compare the coefficient estimates with the corresponding values obtained from wind tunnel tests 
an important distinction between the wind tunnel model and the estimated HARV parameters had to be made. The 
tunnel database models a "basic airframe" normal force coefficient (C NJ J which represents the linear combination 
of the normal force aerodynamic bias (C No ) and the effect of the lift curve slope. Similarly, a basic airframe 
moment coefficient (C m bi ) is presented which involves a linear combination of pitching moment aerodynamic bias 
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{Ca J with the effects of the pitching moment coefficient due to alpha. These coefficient representations used in the 
wind tunnel modeling are shown below: 

Quw = Cl N c * * (2.14) 

= + (2-15) 

All other model parameters were compatible for straight forward comparison with the estimated HARV 
derivatives with the exception of the thrust vectoring control derivatives for which, of course, aerodynamic wind 
tunnel data were not available. Since the estimated coefficients of C No , C Na , C ln() and C [:KI were not directly 
presented in the tunnel database these estimates were combined in the aforementioned manner to present the 
estimated "basic airframe" normal force and moment coefficients for the HARV. 

The estimates of the normal force coefficients are shown in the left side of Figures 5-7. In general, the 
normal force coefficients were determined with consistency at each angle of attack with clear trends over the entire 
alpha range. In addition, the Cramer-Rao bounds closely model the scatter achieved by some of the estimates, 
especially at lower angles of attack. However, many of the bounds show evidence of underestimation. Overall, 
it is clear that each flight maneuver tested is equally capable of concisely estimating most of the normal force 
coefficients. However, there are two main exceptions to this conclusion. Examination of the estimates for C N5pv 
in Figure 7 clearly indicates that the PPSSTV maneuvers, analyzed with all elevator derivatives held constant to 
wind tunnel values, result in more consistent estimates over the entire alpha range. The estimates for the normal 
force coefficient due to pitch rate (C Nq ) in Figure 7 also show that the PPSS and PPSSTV maneuvers yield more 
consistent results with very little scatter over the range of alpha investigated. The SSI MD maneuver results, on 
the other hand, indicate considerable scattering in this estimate This is not surprising since the pilot inputs of the 
PPSS and PPSSTV class maneuvers involve much more sustained excitation in pitch rate than those of the SSI MD 
OBES inputs as shown in Figure 8. 

The results for the estimates of the pitching moment coefficients do not show, in general, the same 
consistency found in the estimates of the normal force coefficients although most estimates still indicate a detectable 
trend. In fact, the results show a severe scatter of the estimates, particularly at high angles of attack. However, 
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closer examination indicates that the majority of the scatter in these coefficients is due to the estimates obtained using 
the SSI MD maneuver data. In fact, it appears that the pitching moment estimates found through analysis of the 
PPSS and PPSSTV data are much more consistent at each angle of attack. This is especially evident in the estimates 
for and C^. Again, this may be due to the highly sustained pitch rate excitation produced by these maneuver 
types. 

The PID estimates for C N b , show good consistency with the wind tunnel data at low angles of attack. At 
higher values of alpha the overall "basic airframe" normal force estimates are lower than the corresponding tunnel 
data. The consistency and low scattering present in the C No and C Na estimates may suggest that the estimates from 
flight data are more realistic than the wind tunnel predicted values. The "basic airframe" pitching moment estimates 
(C ra J also indicate a close match with the tunnel data at lower alphas. However, again there is a substantial scatter 
in the estimates at higher angles of attack. The consistent estimates from the PPSSTV maneuver indicate a higher 
negative basic airframe pitching moment than was predicted by the wind tunnel tests at higher alphas. Comparisons 
between PID estimates and wind tunnel data for the elevator normal force (C N5e ) and moment coefficients (C :;i6 J 
are also shown in Figure 5. Overall, the trends of these estimates match the trend relative to wind tunnel data. 
However, estimates for C N6e at a =40" are consistently of lower magnitude than the wind tunnel data. 

The estimated normal force and moment coefficients for symmetric ailerons (C N5s „ and trailing edge 
flaps (C NSlef , C^f) are compared to wind tunnel data in Figure 6. There is only a limited agreement between the 
PID estimates and the wind tunnel estimates. 

Figure 7 shows instead the estimates of the normal force and pitching moment coefficients associated with 
the pitch vanes of the thrust vectoring system (C N6pv , C rr/1 ..,) as well as the dynamic coefficients pertaining to pitch 
rate (C Nq , C^). Among the longitudinal maneuvers, the PPSSTV maneuvers, performed with OBES off and the 
TVCS on, were particularly effective in extracting estimates for C NSpv and which reveal a clear and precise 
trend with limited scatter as opposed to the other maneuvers analyzed. The trend is made more noticeable by the 
least squares fits shown in the top of Figure 7 which were generated using PPSSTV estimates only. The results for 
C Nq and may again be seen in the bottom of Figure 7. As previously mentioned, the estimates for C Nq and 
resulting from analysis of the PPSS and PPSSTV maneuvers are much more consistent over the entire alpha range 
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tested with very little scatter. Consequently, the least squares fits shown on these plots were generated using the 
PPSS and PPSSTV estimates only. Overall, flight tests indicate a much higher normal force due to pitch rate than 
was predicted during wind tunnel tests. Similarly, the least squares fit for the PPSS and PPSSTV data indicate a 
trend for pitch damping (C^ that is slightly higher in magnitude than that predicted by the wind tunnel data. 
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Chapter 3: Estimation of the Lateral-Directional Aerodynamic Parameters 


3.1 Mathematical Modeling of the Aircraft Dynamics 

Modeling of the aircraft lateral dynamics begins with the polar coordinate form (a,p,V) of the classical 
6 DOF, non-linear equations of motion shown in Refs. [9,11]. The body axis lateral-directional equations of motion 
are given by: 

P=psino-rcosa + -^^Cj. wind +-^ (cosPcos0sin«J>) + 


-2? (-sinP (cos0cos<|>sina-sin6cosa) ) +-^cosasinP 


(3.1) 


Ixfi Ixy&~^xz^ ~ [<?1 (Ty + (q 2 -Z ) Iy z +P<jlxz P-TX^y] g 


(3.2) 


I z i xz& ~ [pq( I x ~Iy) + (P Q ) Ixy*P*-Iyz 9T^X xz ] ^ 


(3.3) 


<J> = p + tan0(rcos<|) + <jsin4>) 


where : 

Cyvin^CyCOsP+CDSinP (3.5) 

C D =C A cos (a) +C w sin(a) (3.6) 

In addition, adjustments were built into pEst to include the effects of angular rates and accelerations on the computed 
data to account for accelerometers which are non-coincident with the aircraft CG. The additional estimated 
parameter is necessary to account for the accelerometer instrumentation bias as shown: 
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a y = C Y + ( x a y ( z a y ~ z cg)i^ ^2 ^a y “/pg) (p 2 +r 2 ) + a yBia ff (3.7) 

Similarly, the flow angle measurements for p recorded by the aircraft instrumentation are also affected by angular 
velocities and sidewash if the flow vane is not coincident with the aircraft CG. Consequently, the computed flow 
angle P is adjusted within pEst as follows: 

P = kiW+(z r z cg )£-(x r x cg )±) (3.8) 

Sensor and CG locations above, behind and to the right of the reference point are considered positive 3 . The 
reference point is located approximately 5 ft in front of the aircraft nose. 

The linear buildup used to model the total, non-dimensional aerodynamic force and moment coefficients 
within pEst are as follows: 

c y“ c y 0 + c y p P + + C Y,^dht + C Yb 6r + 

2^) p * c ^) r <3 - 91 

C, ■ cv C 1,V * c hfi a * c h^ l “ mc * c i.. ir * 

c, /^_\p + c, /— —\r (3.10) 

J ^2W?r M 2VR] 

c 0 = C no - C fl> P + c ^6 a + c^bdht * C^br + 

° n {'TVR) P + C ” r ( 2 V7?) r t 3 - 11 ) 

It should be noted that the static buildup for the rolling and yawing moment coefficients do not include the 
derivatives modeling moment effects due to yaw vane deflections. These moment derivatives are found in the total 
rolling and yawing moment equations: 
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Mx ~ qSbCjR + [C hy byv]R 


( 3 . 12 ) 


M z = qSbC n R + 6q(ERPM(I xg ) ) + [C^Jiyv] R ( 3 . 13 ) 

The second term on the right side of Eq. (3.13) is related to the gyroscopic effects of the engines. The term ERPM 
denotes engine revolutions per minute and I xe represents the moment of inertia of the engines. The factor of 6 is 
simply a conversion factor from revolutions per minute to degrees per second. A standard sign convention for 
control surface deflection was employed: positive rolling control surface deflections induce a positive roll rate, 
positive rudder deflection is defined as trailing edge left. The overall yaw vane deflection was considered positive 
for deflection of the exhaust plume to the right. 


3.2 Estimation Procedures 

Three types of lateral maneuvers were available for PID purposes. The SSI MD maneuvers consisted of 
independent lateral control surface doublets generated by the On Board Excitation System (OBES) for the ailerons, 
differential horizontal tail, rudder, and yaw vanes. The OLAT runs included a series of 5° aileron and 10° rudder 
doublets applied by OBES. Recall that the OBES is instrumental in generating the independent control surface 
deflections necessary for more accurate estimates by minimizing correlated control surface inputs. Each of the 
previous maneuvers were performed by the OBES with the TVCS active. The PYRS maneuvers, which were not 
executed using OBES, were separated into two subclasses pilot yaw/roll steps without thrust vectoring (PYRS) and 
with thrust vectoring (PYRSTV). The PYRS maneuvers consisted of a series of 1.5" rudder steps followed by a 
series of 2.5" lateral stick steps both alternating in sign The PYRSTV inputs consisted of 1.5" rudder steps 
followed by 1.5" lateral stick steps of alternating sign Representative plots of the control surface inputs for the 
maneuver types described above may be seen in Figure 9 recorded at an alpha of 20°. For each maneuver class, 
the lateral inputs were designed to minimize longitudinal dynamic responses. All maneuvers were initiated under 
trimmed, straight and level flight conditions with all flight data sampled with a frequency of 40 Hz. 

Concerns regarding a possible dynamic instability and overestimation in magnitude of estimated derivatives 
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due to correlations in the flight data induced a need for a model simplification. The presence of correlation among 
static and/or dynamic flight data implies that a dependency exists among them often leading to difficulties in 
extracting model parameters that describe their separate effects on the overall system. These correlations often exist 
among control surface deflections which are dependently deployed with respect to one another. An analysis of the 
HARV lateral flight data for the PYRS maneuver clearly indicates a high correlation between aileron and differential 
horizontal tail inputs when the standard flight control system is active. There is also a dependency between roll and 
yaw control inputs applied by the aileron-rudder interconnect (ARI) of the standard F/A-18 flight control system. 
These highly correlated inputs are clearly indicated in Figure 10. As a result, all differential horizontal tail 
derivatives (C Y5dhl , C l5dhl and C n5iht ) were held constant at wind tunnel values during the analysis of the PYRS data 
so that more accurate aileron coefficients could be obtained. The operation of the ARI during these pilot yaw/roll 
step maneuvers threatened to introduce difficulties in extracting the cross-coupling control derivatives C Y6s , C n6j 
and C 16r . Analysis of the PYRS maneuvers at 20" angle of attack achieved best results when these cross-coupling 
control derivatives were held constant at the estimates from wind tunnel analysis. However, PYRS data analysis at 
30” angle of attack and higher required that only the cross-coupling control derivative C n6l be held constant at 
appropriate wind tunnel values. 

Similarly, examination of the PYRSTV flight data also indicates a high correlation between aileron and 
differential horizontal tail inputs as well as a dependency between roll and yaw control inputs applied by the ARI 
of the standard F/A-18 flight control system. In addition, when the TVCS is active there is an additional correlation 
between rudder and yaw vane inputs. A sample time history showing these correlations may be seen in Figure 1 1 
Best results were obtained during the analysis of the PYRSTV flight data by holding all differential horizontal tail 
derivatives (C YMtl , C 15dhl and C nMhl ) as well as all rudder coefficients (C Y5r , C 1Sr and C nSr ) constant at appropriate 
wind tunnel values in order to obtain better estimates for the aileron and yaw vane derivatives. All wind tunnel 
values were based upon angle of attack and Mach for the maneuvers. 

All mass properties (x,,, y cg , z C8 , I x , I y , I 2 , I u and mass) were taken to be unique values for each individual 
maneuver. These values were obtained from flight data with an average value over the appropriate maneuver time 
being used in the mathematical model found in pEst. Finally, all longitudinal states and responses found in the 
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aircraft equations of motion, including the aircraft flight velocity V, were replaced with measured flight data 
throughout the estimation process. 

The following static and dynamic lateral-directional derivatives were estimated for the NASA F/A-18 
HARV: C Yo , C Y5a , C Y6 dhi, C Yfir , C Y6yv , C Ypf C Yr , C ]0 , C l{3 , C^, C 16dht , C 15r , C [5yv , C ]p , C lr , C Q0 , C n p, 

C n5l , C nMht , C n5r , C n6yv , C np , C.. The additional parameter aye,, was estimated in order to model the 
measurement bias for the lateral accelerometer aboard the HARV. A complete listing of which aerodynamic 
coefficients were estimated, as well as which derivatives were held constant at wind tunnel values, for each class 
of maneuver may be seen in Table 3. 

In terms of the lateral-directional cost function (J), the objective was to minimize the difference between 
the measured and computed values of p, p, r, <\> and ay as shown below. 




VL y ' 
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A 0 0 0 0" 


-P(t k ) 


0 B 0 0 0 




0 0 C 0 0 




0 0 0 D 0 




0 0 0 0 E 


~a y (t k ) 


L 


p(t k ) -p(t k ) 
r(t k ) -f( t k ) 
<t> ( t k ) ( c k ) 


(3.14) 


Here the terms (3, p, r, 6 and a, are the responses measured from the aircraft while the corresponding responses 
computed in pEst are denoted by (3, p, f, 4> and a,. The diagonal matrix within the cost function represents the user 
defined weighting matrix (W,). A trial and error approach was used to determine the best set of weightings for each 
class of maneuver investigated. A listing of the A, B, C. D and E elements of the weighting matrix, as was used 
for each class of maneuver, may be seen in Table 4 


3.3 Results of the Parameter Estimation Process 

The results for the lateral-directional stability derivatives of the NASA F/A-18 HARV are presented in 
terms of time history comparisons between measured and computed aircraft responses as well as estimates for the 
aerodynamic coefficients which are shown in companson with wind tunnel data. Overall, estimates were obtained 
for the aerodynamic subcomponents of the total roll, yaw and side force coefficients (C,, C D and C Y ). 
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Comparison between computed aircraft responses, generated using the estimated aerodynamic model, and 
measured time histories indicate that an accurate representation of the HARV lateral-directional flight data responses 
have been obtained even at higher alphas. Figure 12 offers a representative sample of the most important time 
history fits of (3, p, r and 4> at an alpha of 30° for the SSI MD class maneuver. 

Estimation results for the main lateral-directional stability derivatives are presented in Figures 13-16 over 
a flight regime of 5°- 60° angle of attack. Each estimate is shown with its corresponding Cramer-Rao bound 
represented graphically as a vertical bar. Recall that the Cramer-Rao bound is a computed lower bound of the 
standard deviation of each estimate which, under ideal estimation conditions, should accurately model the scatter 
of the estimates at each angle of attack 6 . Some considerations are needed for a better understanding of results in 
terms of the Cramer-Rao bounds. In fact, recall that the determination of the Cramer-Rao bounds assumes unbiased 
estimates with system and measurement noise being described as a Gaussian, white sequence. Considering a case 
where these assumptions are truly valid the time history residuals would also be white, random sequences resulting 
in a set of Cramer-Rao bounds accurately modeling the scatter of the estimates. However, such ideal test conditions 
are only guaranteed if flight data is obtained through computer simulation. As expressed in Ref. [6], analysis of 
real flight data may result in biased estimates often caused by modeling inadequacies or poor system excitation 
Furthermore, colored system or measurement noise will also affect the computed residuals. As a result, in the case 
of residuals correlated in this manner, optimistically underestimated Cramer-Rao bounds are often calculated causing 
an inadequate model for the scatter of the estimates. 

A variety of techniques to "correct" the Cramer-Rao bounds are discussed in Ref. [6]. One such technique 
involves multiplying the bounds by a factor, often 5-10, to compensate for the effects of the colored residuals. A 
more precise method discussed in Ref. [6] entails determining a unique correction multiplier for each bound through 
examination of the power spectral density for the response residuals generated by the estimates. Note that the 
Cramer-Rao bounds presented herein are not adjusted. 

Overall, estimates of the static derivatives show good consistency with well defined trends. In addition, 
some of the Cramer-Rao bounds are successful in modeling the estimate scatter at the lower angles of attack 
However, the majority of the bounds are clearly underestimated. 
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An interesting trend was determined, as seen in Figure 13, for the dihedral effect (C ip ) in which the 
estimates start at a low magnitude at an angle of attack of approximately 5 -10* and clearly increase to a maximum 
magnitude at about 20* alpha. The magnitude of the dihedral effect decreases steadily thereafter. The results for 
the aircraft yawing coefficient due to sideslip (C nP ) indicate that the HARV is directionally stable in yaw with C„ P 
positive up to approximately 30’ angle of attack. However, estimates indicate that the aircraft becomes unstable 
in yaw at alpha > 30’, as is indicated with C op being negative. 

Figures 13-14 indicate an expected degradation in the rolling effectiveness of the ailerons (C,J and 
differential stabilizer (C 1M1U ) as the angle of attack increases. This may be due to slow speeds and flow separation 
which occur under high alpha flight conditions. Similarly, as seen in Figure 14, the yaw control of the rudder 
(0,^) was found to become less efficient at higher angles of attack. This could be caused, most likely, by the 
rudder being in the wake of a highly turbulent, separated airflow from the wing under such flight conditions. 
Although the estimates for the cross-coupling control derivatives (C^, C n5dhI and C 16r ) do not show the same 
consistency as the primary control derivatives (C 15 ,, C 15dhl and C QSr ) a clear trend is still indicated. Finally, as it 
can be seen in Figure 16, the yaw control effectiveness of the thrust vectoring vanes (C^) aboard the HARV 
remains almost constant over the entire flight regime with only a slight variation in magnitude. 

The main lateral-directional side force coefficient estimates with the most defined trends may also be seen 
in Figure 16. Again, the majority of the scatter is obtained at higher alphas with the estimates for the SSI MD 
maneuvers. As expected, the side force effectiveness of the yaw vanes (C YSyv ) remains nearly constant over the 
entire flight regime investigated. The results for the side force effectiveness of the rudder (C Y5r ) indicate a loss 
of rudder effectiveness as alpha increases. In addition, the scatter in the estimates for C Y5r at higher alphas indicate 
a possible rudder control reversal under these flight conditions. However, such a control reversal is clearly not 
present in the rudder yaw control results for the C n5r estimates previously examined in Figure 14. 

Among the dynamic coefficients estimated, the most consistent results were obtained for C lp and as can 
be seen in Figure 15. Overall, these coefficients indicate a damping effect against roll and yaw respectively. The 
majority of the scattered dynamic derivative estimates are represented by results obtained using the OLAT inputs. 
This may be due to a lack of proper excitation in roll and yaw resulting from the OLAT maneuvers in the estimation 
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of these coefficients. In addition, some scatter of the SSI MD estimates is present at higher alphas. 

A comparison between the main lateral-directional aerodynamic estimates from flight data with the 
corresponding values in the NASA F/A-18 HARV wind tunnel database was also conducted. These final 
comparisons may also be seen in Figures 13-16. Each plot contains a least squares curve fit which has been 
generated through the estimates for comparison with the appropriate wind tunnel values when available. Note that 
the F/A-18 wind tunnel database does not include values that are directly comparable to C ip , C nP and C Yp . In fact, 
the wind tunnel data linearly combines the effects of each sideslip coefficient with the appropriate roll, yaw or 
sideforce aerodynamic bias coefficient into one overall term. In addition, no compatible wind tunnel data was 
available for direct comparison with the yaw vane coefficient estimates. 

Figure 13 shows a comparison between the estimates, and corresponding tunnel data, for the rolling and 
yawing coefficients due to aileron deflection (C 15l and C n5a ). The estimates for C l6 , are in close accordance with 
the listed tunnel data whereas the results for C n5 , are consistently lower in magnitude than the tunnel data. This 
trend indicates less adverse yaw than expected due to aileron deflection at higher angles of attack. 

Comparisons for the rudder roll and yaw estimates are presented in Figure 14. The estimates for C n5r 
match the tunnel data with good accuracy. However, estimates for adverse roll due to rudder deflection (C l5r ) are 
close to wind tunnel estimates at the low and high ends of the flight regime between 10°- 60° alpha. In the region 
of 30°- 40° alpha the estimates for C 15r are consistently higher than those available in the wind tunnel database. 

Estimates for C 15dhI , as seen in Figure 14, are consistently higher than those of the wind tunnel data. 
Similarly, resulting estimates for C nSdht show that, as angle of attack increases, the flight data indicate a higher 
magnitude of adverse yaw due to differentially deflected stabilizers than was previously modeled in the wind tunnel 
database. 

The least squares fit of the dynamic coefficient estimates of C, p and C m are shown in Figure 15. Both fits 
show good agreement with the wind tunnel data. Similar results were obtained for the least squares fits of the 
dynamic coefficient estimates of and C np . 

Comparisons of the estimates for C Y * with wind tunnel data may be seen in Figure 16. Again, the least 
squares fit through the estimates agrees with the wind tunnel values, especially at low angles of attack. However, 
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closer examination shows that for oc^40° a major divergence from the wind tunnel data is encountered in which 
the flight data indicates a possible rudder control reversal. Again, this control reversal was not indicated by the 
C n6r estimates. 
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Chapter 4: Conclusions 


This paper has presented the longitudinal and lateral-directional stability and control derivatives 
obtained for the NASA F/A-18 HARV from flight data using the Maximum Likelihood estimation technique 
with a Newton-Raphson minimization scheme over a flight regime of 5°- 60' angle of attack. The 
mathematical model used in this investigation included classical static and dynamic derivatives estimated through 
the analysis of a variety of PID maneuvers. 

Resulting time history comparisons, along with the examination of the derivative estimates, have shown 
that the final estimates are adequate representations of the dynamic system. The estimates of the normal force 
coefficients, in general, indicated clear trends, exhibited very good consistency along with a limited level of 
scatter, and a reasonable size for the Cramer-Rao bounds. The estimates of the pitching moment coefficients 
were found to be best extracted from flight data pertaining to the PPSS and PPSSTV maneuvers with their 
inherent ability to sustain excitation in pitch rate. Consequently, these maneuvers also produced the most 
consistent results for the dynamic coefficients C Nq and C^, in comparison with estimates obtained from other 
maneuvers. The normal force and moment coefficients due to pitch vane deflection were found to be best 
identified using the PPSSTV class maneuvers. Overall, very good agreement of the estimates with the tunnel 
data was obtained at lower angles of attack with some deviations apparent at higher alphas. The majority of the 
rolling and yawing moment coefficients, indicated good consistency with well defined trends and limited scatter 
present mostly at higher alphas. However, although discernible trends are still evident, less consistency in the 
estimates was achieved for the cross-coupling control derivatives C^, C nf<Ull and C 15r than the primary control 
derivatives and C R{r . 
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* For the SSI MD maneuver at an alpha of 60° E = 250 
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Table 3: Breakdown of aerodynamic coefficients estimated for each class of maneuver analyzed. 
(✓ = Parameter was estimated, W/T = Coefficient held constant at wind tunnel value, 
N/A = Surface was not active and thus no estimate was obtained) 



* Indicates parameter was estimated for all alpha higher than 20* 
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